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We study the effect of laser phase noise on the generation of stationary entanglement between 
an intracavity optical mode and a mechanical resonator in a generic cavity optomechanical system. 
We show that one can realize robust stationary optomechanical entanglement even in the presence 
of non-negligible laser phase noise. We also show that the explicit form of the laser phase noise 
spectrum is relevant, and discuss its effect on both optomechanical entanglement and ground state 
cooling of the mechanical resonator. 
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I. INTRODUCTION 

The study of optomechanics of micro- and nano- 
cavities has recently sparkled the interest of a broad sci- 
entific community due to its different applications, rang- 
ing from sensing of masses, forces and displacements at 
the ultimate quantum limits [TJ [2] , to the realization of 
quantum interfaces for quantum information networks 
[3]-[6], up to tests of the validity of quantum mechan- 
ics at macroscopic level [TJ [8] . The possibility to detect 
genuine quantum behavior in cavities characterized by 
an appreciable radiation pressure interaction between a 
light mode and a mechanical resonator was first pointed 
out by Braginski and coworkers [9] in the context of the 
interferometric detection of gravitational waves. In the 
last years however many different schemes have been pro- 
posed for the detection of quantum mechanical effects in 
such systems, such as continuous variable (CV) entangle- 
ment between cavity modes and/or mechanical modes, 
squeezed states of the light or the mechanical modes, 
and ground state cooling of the mechanical modes (see 
Ref . [10] for a review) . These schemes involve cavities and 
resonators at the micro- or nano-level rather than the 
macroscopic scale of gravitational wave detectors, and 
profit from the tremendous progress in micro- and nano- 
fabrication techniques which has provided novel oppor- 
tunities to engineer optomechanical devices. Some exam- 
ples are toroidal optical microresonators [2 , Fabry-Perot 
cavities with a movable micromirror [HJ[12], a semitrans- 
parent membrane within a standard Fabry-Perot cavity 
[T3HT6] . suspended silicon photonic waveguides [T7UT9] . 
SiN nanowires evanescently coupled to a microtoroidal 
resonator [20], adjacent photonic crystal wires [21], nano- 
electromechanical systems formed by a microwave cavity 
capacitively coupled to a nanoresonator [22-24 , atomic 
ensembles interacting with the mode of an optical cavity 
containing it [25H27] . 

Here we focus on the possibility to generate robust CV 
entanglement between an intracavity optical mode and 
a mechanical resonator mode in the stationary state of 



the system, which has been predicted in [28 (see also 
Refs. [29H3T] where the problem has been revisited, and 
Ref. [32 for its extension to the microwave cavity case). 
This steady-state entanglement could be very useful for 
quantum communication applications because it is end- 
less, and robust against thermal noise. The above analy- 
ses however have not taken into account a possible techni- 
cal limitation, associated with the fact that both the am- 
plitude and especially the phase of the laser driving the 
cavity are noisy quantities. Laser phase noise could be 
in fact very dangerous for optomechanical entanglement, 
which is just the existence of strong quantum correla- 
tions between the fluctuations of cavity field quadratures 
at an appropriate phase and the position and momentum 
fluctuations of the mechanical resonator. 

The effect of laser phase noise on ground state cool- 
ing of the mechanical resonator has been already dis- 
cussed in Refs. [33-36 . Ref. [33] showed that phase noise 
acts on the mechanical resonator as an additional heat- 
ing noise proportional to the intracavity field amplitude, 
that may represent a serious obstacle for ground state 
cooling. Ref. [35 showed that what is relevant both for 
cooling and also for coherent state transfer between op- 
tical and mechanical modes is just the phase noise spec- 
trum at the mechanical resonance frequency. If such a 
noise value is not too large, cooling is still possible as, in 
fact, has been confirmed experimentally in Refs. [37] [38] 
which approached, and Ref. [39] just reached, the ground 
state limit. Therefore the explicit form of the laser phase 
noise spectrum is relevant, and simply assuming a white 
phase noise tends to overestimates its effect because such 
a noise is strongly colored and decays significantly at 
MHz frequencies. 

Here we extend the analysis to stationary optomechan- 
ical entanglement, in order to establish to what extent 
laser phase noise may affect its realization. Our analy- 
sis is based on a quantum Langevin equation treatment 
and generalizes previous approaches in various aspects. 
In particular it can deal with various examples of phase 
noise spectra, even though we shall focus on a bandpass 
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filter form of the laser noise spectrum. This is in fact 
the typical case in current experiments because phase 
noise decays to zero at large frequencies, and it is negli- 
gible at low frequencies due to laser-cavity locking. We 
show that laser phase noise has an appreciable effects 
on the achievable entanglement, but nonetheless signifi- 
cant stationary entanglement is still achievable employ- 
ing currently available stabilized lasers. We also derive 
approximate analytical expression for the stationary op- 
tomechanical entanglement illustrating how laser phase 
noise affects its experimental realization. We shall also 
briefly reconsider the problem of mechanical ground state 
cooling. 

The paper is organized as follows. In Sec. II we present 
the model Hamiltonian and we show how one has to mod- 
ify the standard quantum Langevin treatment in order 
to include in a non-perturbative way the effects of laser 
noise. In Sec. Ill we study the linearized dynamics of 
quantum fluctuations around the appropriate classical 
steady states. In Sec. IV we show the results for sta- 
tionary optomechanical entanglement and cooling to the 
mechanical ground state. Sec. V is for concluding re- 
marks. 



II. MODEL 

We consider a generic cavity optomechanical system 
in which a mechanical resonator with frequency cj m is 
subject to a force proportional to the photon number 
of an optical cavity mode with frequency Co> c , which is 
driven by an intense, but noisy laser. The corresponding 
Hamiltonian can be written as [2j [TUJ [40j [41] 

H = huj c a^a + -huu m (p 2 + q 2 ) — HGqo) aq 

+ ih£(t)(aU- { ^ ot+ *W - ^ot+Ht)]^ (i) 

The first term describes the energy of the cavity mode, 
with annihilation operator a ([a,af] = 1), while the sec- 
ond term gives the energy of the mechanical resonator, 
described by dimensionless position and momentum op- 
erators q and satisfying the commutation relation 
[QiP] — i- The third term is the optomechanical in- 
teraction, with single photon optomechanical coupling 
strength 

Go = -(p) J^, (2) 
\ ax J V Triu m 

where (duj c /dx) is the change in cavity frequency per dis- 
placement and m is the effective mass of the mechanical 
mode [42] . 

The last term describes the cavity driving by a laser 
which is generally assumed to possess both phase and in- 
tensity fluctuations. The parameter uq denotes the laser 
average frequency and <p(t) is the zero- mean fluctuating 
phase, while £{t) = So + e(t) is related to the laser am- 
plitude and e(t) describes the real, zero- mean amplitude 



fluctuations of the laser. The statistical properties of <j>(t) 
and e(t) will be specified later on, while the mean ampli- 
tude is given by £q = ^JlnP '/fouio, where P is the input 
laser power and k is the cavity loss rate through its input 
port. 

The Hamiltonian of Eq. describes a wide variety of 
cavity optomechanical systems, with different cavity ge- 
ometries and mechanical elements. In systems as Fabry- 
Perot cavities with a movable micro-mirror [TTJ [12j [431 
|45] , or with a semitransparent membrane inside [T3l - 
[16], or in radially vibrating toroidal microcavities [46] . 
optomechanical coupling is provided by radiation pres- 
sure. In other optomechanical devices coupling is in- 
stead provided by the transverse gradient force, such as 
in suspended silicon photonic waveguides [T7HT9] . SiN 
nanowire evanescently coupled to a microtoroidal res- 
onator [20], and in "zipper" cavities formed by two adja- 
cent photonic crystal wires [21]. The same Hamiltonian 
applies also to nanoelectromechanical systems formed by 
a microwave cavity capacitively coupled to a nanores- 
onator, such as in [22-24 , and in such a case the noisy 
laser describes the phase and intensity fluctuations of the 
microwave driving source. Finally Eq. also applies to 
systems where a mechanical collective degree of freedom 
of an atomic ensemble interacts with an optical cavity 
containing it [25-22]. In all these devices one always has 
many cavity and mechanical modes, but one can adopt 
the Hamiltonian of Eq. ([I]) whenever one can restrict 
to single cavity and mechanical modes. This is justi- 
fied when the cavity free spectral range is much larger 
than the mechanical frequency u m (i.e., not too large 
cavities): in such a case the input laser drives only one 
cavity mode and scattering of photons from the driven 
mode into other cavity modes is negligible [40 . One can 
restrict to a single mechanical mode when the detection 
bandwidth is chosen so that it includes only a single, iso- 
lated, mechanical resonance and mode-mode coupling is 
negligible [47] . 

Quantum Langevin equations 

For a full description of the system dynamics it is nec- 
essary to include the fluctuation-dissipation processes af- 
fecting both the optical and the mechanical mode. They 
can be taken into account in a fully consistent way [41], 
and one gets 

q = u m p, (3a) 

p -Umq-lmP + Goata + Z, (3b) 
a = — [k + \uo c ]a + iG$aq 

+£(t)e- i l Uot+ *M + V^a in , (3c) 

where ai n is the vacuum input noise, whose nonzero cor- 
relations are given by [48] 

(a in (t)al(t')) = [N(oj c ) + l]5(t-t'), (4a) 
(al(t)a in (t')) = N(uj c )S(t-t'), (4b) 
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FIG. 1: Scheme of the cavity optomechanical system under 
study monitored by an homodyne detection apparatus. 



a consequence Eqs. ([3| become 

q u m p, (6a) 
p = -uj m q - 7 m p + God a + f, (6b) 
a = —Ka — i(AQ — <j) — Goq)a + £(t) + V^a[ n (t)(fic) 

where Ao = uo c — ujq is the detuning of the cavity 
mode from the average laser frequency, and a- m (t) = 

/,\ ioJo(t— to)+i f* dt'Mt') • ,.11 • , 

ai n (t)e Jt o ^ v is still vacuum input noise, 

possessing the same correlation functions of a[ n (t) [see 
Eqs. Q]. Therefore in the frame rotating at the fluc- 
tuating frequency, laser amplitude noise acts as additive 
noise on the cavity modes, while laser frequency noise is a 
multiplicative noise, affecting the cavity field in the same 
manner of the fluctuations of the resonator position q. 



III. LINEARIZED DYNAMICS OF THE 
FLUCTUATIONS 



with N(ou c ) = (exp{/kj c /k^T} — 1) _1 the equilibrium 
mean thermal photon number (&b is the Boltzmann con- 
stant and T is the temperature of the reservoir). At opti- 
cal frequencies fyjj c /k^T ^> 1 and therefore N(uj c ) « 0, so 
that only the correlation function of Eq. (|4^l) is relevant. 
We have assumed for simplicity the ideal case of a single- 
ended cavity, so that the total cavity decay rate coincides 
with the loss rate through the input port n. The mechan- 
ical mode is affected by a viscous force with damping rate 
7 m and by a Brownian stochastic force with zero mean 
value £(£), obeying the correlation function [48| [49] 

«<«<'» -£/£^»H5Sr)4 

(5) 

The Brownian noise is a Gaussian quantum stochas- 
tic process and its non-Markovian nature (neither its cor- 
relation function nor its commutator are proportional 
to a Dirac delta function) guarantees that the QLE of 
Eqs. ^ preserve the correct commutation relations be- 
tween operators during the time evolution [4T] . 

We are interested in the entanglement properties of 
the stationary state of the system, which are determined 
by the quantum correlations between the mechanical and 
optical field quadratures. These quadratures can be de- 
tected by homodyning the cavity output and an addi- 
tional weak field probing the mechanical element (see for 
example Ref. [28 ). The local oscillator for the homodyne 
detector is provided just by the noisy driving laser (see 
Fig. [TJ and this means that all detected quantities are 
referred to the frame rotating at the fluctuating instan- 
taneous frequency uo + <j)(i). Passing to this randomly 
rotating frame, the cavity field operator transforms ac- 
cording to a(t) — » a(t) exp j— \ujo(t — to) — i f* Q dt' 4>(t')^ , 

where to — » — oo is the time instant at which we fix the 
phase reference for the cavity field by taking £{t) real; as 



Achieving stationary optomechanical entanglement 
means establishing strong quantum correlations between 
the steady-state fluctuations of the position and momen- 
tum of the resonator, and the intracavity field quadrature 
fluctuations. As shown in Refs. [28, 29 , this is attained 
when the effective coupling between these fluctuations is 
strong, which is realized when the intracavity field is very 
intense, i.e., for high-finesse cavities and enough driving 
power. Therefore we focus onto the dynamics of the fluc- 
tuations around the classical steady state. 

When the system is stable, it reaches a steady state 
which, in the absence of laser noise, is characterized by 
the cavity mode in a coherent state, and the mechanical 
resonator at an equilibrium position shifted by a quantity 
proportional to the stationary intracavity photon num- 
ber. One may expect that, due to laser phase noise, this 
classical steady state is modified: in particular one ex- 
pects that the phase of the intracavity coherent state 
slowly becomes completely random. This is the classical 
steady state assumed in Refs. [35, 36 , i.e., a cavity field in 
a coherent state with completely random phase but with 
time-independent photon number, so that the shift of 
the equilibrium position of the resonator is not changed. 
Our treatment adopts however the frame rotating at the 
fluctuating instantaneous laser frequency uo + </>(£), dif- 
ferently from Refs. [35] [36] which adopt the frame rotat- 
ing at ooo- In such a fluctuating frame the phase of the 
classical stationary coherent state is not random, and its 
amplitude a s is given by the standard expression valid in 
the absence of laser noise, a s = £o/(k + iA). This latter 
equation is actually an implicit nonlinear equation for a s 
because 

A = Ao-Go 9s = Ao-^^, (7) 

the effective cavity detuning, depends upon |a s | 2 . As- 
suming such a time-independent classical steady state 
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is equivalent to assume that both phase and amplitude 
noise of the driving laser affect only the quantum fluc- 
tuations of the system and not its classical stationary 
state. In fact, inserting such a steady state solution into 
Eqs. (J6|, the exact (nonlinear) QLE for the fluctuations 
become 

Sq u m Sp, (8a) 
Sp = —uj m Sq — lm.Sp + Go(a s 5a) + a* 8a) + Sa^Sa -(8&) 
Sa = —[K-\-iA]5a + iGoa s 5q + i[Go5q + <j)]5a 

+i^a s + e + Vatican- (8c) 

As discussed above, robust optomechanical entanglement 
can be generated when the effective coupling between the 
fluctuations Goa s is large enough, which is best achieved 
when \a s \ ^> 1, i.e., we have a large number of intracavity 
photons. In such a case the system dynamics is well 
described by linearizing Eqs. (J8|, i.e., by neglecting the 
term 8a) 8a in Eq. (J^Jd) and the terms i[GoSq + <j)]8a in 
Eq. so that 

Sq = u m 8p, (9a) 
5p = -cu m 8q - ^ m Sp + G (a s 8a^ + a*8a) + £, (9b) 
8a = — [k + iA]5a + iGoa s &7 + i<Ms + £ + V^a^c) 

Notice that in this way we neglect also the multiplica- 
tive noise term \(j)8a together with the usual nonlinear 
terms, but this is reasonable because when \a s \ ^> 1 
such a term has a negligible effect compared to that 
of the i(j)a s term. It is convenient to rewrite the lin- 
earized QLE in terms of the field quadrature fluctu- 
ations 8X A = (5ae ie *+5ate- ie *) /y/2 and 8Y A = 



Sq 
Sp 

SX A 
SY A 



IkY, 



A j 



(10a) 
(10b) 
2^X^0c) 

(lOd) 



u m Sp, 

-u m 8q - >y m 8p + G8X A + f , 
-nSX A + A5Y A + V2 cos 9 A s 
-k8Y a - A8Y A + G£g 
W2|a s |0 + V^2 sin 6> A £ 

where A = arctan[A/ft], G = GoV^l^sl is the ef- 
fective optomechanical coupling, and we have intro- 
duced the corresponding Hermitian input noise opera- 
tors X% = (a in e WA +a} n e- w *)/y/2 and Y£ = (a in e^ A - 

~a\ n e-^)/iV2. 

Finally we have to specify the statistical properties of 
laser phase and amplitude noise. In currently available 
stabilized lasers, amplitude noise e(t) is negligible with 
respect to phase noise <p(t) and therefore we shall ne- 
glect it from now on, as assumed also in Refs. [33H36] . 
Phase noise instead is typically non-negligible and it is 
responsible for the nonzero laser linewidth Ti. In fact, 
the laser spectrum is given by the Fourier transform of 
the stationary correlation function of the field, that is, 



5l(w) = j dre^ T C(T) = J 



dTe iuJT (exp{i<j>(t+T)-i<j>(t)}}; 

(11) 



(j){t) is well described by a zero-mean stationary Gaussian 
stochastic process, and therefore one can write 



C(t) = ^expjijTd^Oo}' 

= expj-i j\s J\s'(HsW))Y (12) 

If one takes a delta correlated frequency noise 
(<j)(s)<j)(s f )) = 2Ti8(s—s'\ i.e., a flat frequency noise spec- 
trum S^u) = 2Ti, Eq. (121 yields G(r) = e^M, which 
corresponds to a Lorenzian laser spectrum with linewidth 
T^. However, as already pointed out in Ref. [35 , tak- 
ing a flat S^(u) tends to overestimate the effect of laser 
phase noise; in practice the frequency noise spectrum has 
a bandpass filter form and therefore the laser spectrum 
is no more a perfect Lorenzian, but has faster decaying 
tails. A more realistic description is obtained by taking 
the following bandpass filter form of the frequency noise 
spectrum 



n 4 



(n 2 



,212 



2^2 ' 



(13) 



with ft denoting the band center and 7 the bandwidth 
of the frequency noise spectrum, while the noise strength 
still describes the laser linewidth. A flat frequency 
noise spectrum is recovered in the limit Q — >• 00, 7 — >• 00. 
It is straightforward to verify that the frequency noise 
spectrum of Eq. ( 13 ) is reproduced by assuming that the 
frequency noise variable ip = (j) satisfies the following pair 
of Langevin equations 

ip = ne, (i4a) 

6 = -fill;- 70 + y/2Tie(t), (14b) 

where e(t) is a white noise with correlation function 
(e(t)e(t / )) c i = 8(t — t'). If we now attach these two lat- 
ter equations to the linearized QLE of Eqs. ( |l0| ), i.e., we 
treat the variables i/j and as two additional dynamical 
variables of the system, we get that the system dynamics 
in the presence of laser phase noise is fully described by 
the following set of equations in matrix form 



u(t) = Au(t) +n(t), 



(15) 



where u(t) = (8q(t),8p(t),8X A (t),8Y A (t)^(t),0(t)) T is 
the vector of CV fluctuation operators, and n(t) = 



(0, £(£), y/2KX^(i), V2^ n (t), 0, QV2Tie(t)) 1 is the cor- 
responding vector of noises. Moreover, the drift matrix 
A is the 6x6 matrix 



(16) 
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A. Stationary quantum fluctuations 

We are interested in the stationary properties of the 
system: in particular we want to check to what extent 
laser phase noise hinders achieving a steady state with 
distinct quantum properties, in particular characterized 
by CV entanglement between the cavity mode and the 
mechanical element. The realization of such a stationary 
optomechanical entanglement is of particular interest for 
quantum information applications, because it would rep- 
resent a very robust source of persistent entanglement. 

The steady state associated with Eq. (15) is reached 



Eq. (21 ) is a linear equation for V and can be straightfor- 



when the system is stable, which occurs if and only if all 
the eigenvalues of A have negative real part. These sta- 
bility conditions can be obtained for example by using 
the Routh-Hurwitz criteria, and it is possible to verify 
that they are not modified by the presence of laser phase 
noise. Therefore they coincide with those discussed in 
Refs. [28] |29] [50]; in this paper we shall restrict to the 
situation with A > 0, i.e., with a red-detuned laser, and 
in this parameter region the only non-trivial stability con- 



dition is G 2 < (A 2 



a/A. 



The steady state is a zero-mean Gaussian state due to 
the fact that the dynamics of the fluctuations is linearized 
and all noises are Gaussian; as a consequence, it is fully 
characterized by the 6x6 stationary correlation matrix 
(CM) V, with matrix elements 

Vi . = (^(oo)uj(oo) + ^(oo)^(oo)) ^ 



The formal solution of Eq. ( 15 ) yields [29] 

Vij= ds ds'MikWMjtWDuis-s'), (18) 
Jo Jo 

where M(t) = exp(At) and D(s - s') is the diffu- 
sion matrix, the matrix of noise correlations, defined as 
D M (s - s f ) = (n k (s)m(s') + m(s')n k (s))/2. The Brow- 
nian noise £(£) is in general a non-Markovian Gaussian 
noise (see Eq. (J5|), but in the limit of large mechanical 
quality factor Q m — co> m /7 m ^> 1, becomes with a good 
approximation Markovian, with symmetrized correlation 
function 



: <y m (2n + l)5(t-lf), (19) 



where n = (ex.p{huu m /kBT} — 1) _1 is the mean thermal 
phonon number at T. Therefore, the diffusion matrix be- 
comes D(s-s') — DS(s-s'), where D = diag[0, 7 m (2n + 
1), k, k, 0, 2T/0 2 ], so that Eq. (18) simplifies to 



poo 

V= dsM(s)DM(s) T , 
Jo 



(20) 



which, when the stability conditions are satisfied 
(M(oo) = 0), is equivalent to the following Lyapunov 
equation for V 



AV + VA T = -D. 



wardly solved, but its explicit expression is cumbersome 
and will not be reported here. 



IV. RESULTS 



From the solution of Eq. (21) for the stationary CM 
V, we can determine all the quantum properties of the 
stationary state of the cavity optomechanical system. We 
determine in particular the effect of laser phase noise on 
the possibility to achieve optomechanical entanglement 
and ground state cooling of the mechanical resonator. 



A. Entanglement 

The auxiliary variables i/j and do not refer to the 
optomechanical system of interest and therefore we are 
concerned with the reduced 4x4 CM extracted from 



Eq. (20) by neglecting the last two rows and columns. 



This reduced correlation matrix can be expressed in the 
following form 



V 



V A V c 
V B 



(22) 



where Va, Vb and Vq are 2x2 matrices, with Va associ- 
ated to the mechanical resonator, Vb to the cavity mode, 
and Vq describing the optomechanical correlations. A 
convenient measure for CV entanglement is the logarith- 
mic negativity [5lJ [52] , given by 

(23) 



En = max(0, — ln277 



where r] is symplectic eigenvalue of the bipartite system 
and it is given by the equation 



V 



1 

71 



[S(V) - ^(V) 2 -4det V] 



1/2 



(24) 



(21) 



with Y,(V) = det V A + det V B - 2 det V C - 

We now study the behavior of En when the phase noise 
parameters, i.e., the laser linewidth Ti and the center of 
the phase noise spectrum are varied. For simplicity we 
shall always consider a bandpass filter form of the phase 
noise spectrum, i.e., Eq. (13) with 7 = fl/2. We consider 
typical parameter values, i.e., a mechanical resonator 
with cj m /27r = 10 MHz, quality factor Q m = 2 x 10 6 , 
and at T = 0.4 K. We also assume a single photon op- 
tomechanical coupling strength Go = 1 KHz, which is 
achieved for example in a Fabry- Perot cavity with length 
L = 1 mm and with an oscillating micromirror of ef- 
fective mass m ~ 10 ng. However, comparable or even 
larger values of Go are currently achieved in other cavity 
optomechanics setups. 

Typical values of the linewidth of stabilized lasers are 
around F//27T ~ 1 KHz and therefore it is important to 
see if the stationary optomechanical entanglement pre- 
dicted in Refs. [28, 29 is robust against laser phase noise. 
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FIG. 2: (Color online) Contour plot of En versus the input power P and normalized cavity detuning A/cj n 
phase noise), Ti/2n = 0.1 KHz (b), and Ti /2tv = 1 KHz (c). The cavity bandwidth has been fixed at ft/cj m 
also fixed the center of the phase noise band Q/2tt = 50 KHz. See text for the other parameter values. 



for Ti = (a) (no 
= 0.5 and we have 




FIG. 3: (Color online) Contour plot of En versus the input power P and normalized cavity decay rate n/oj m for Yi = (a), 
Ti/2ty = 0.1 KHz (b), and Ti/2n = 1 KHz (c). The cavity detuning has been fixed at A/cj m = 1 and we have fixed Q/2ty = 50 
KHz. See text for the other parameter values. 



In Fig. [2] we show En versus the input power P and nor- 
malized cavity detuning A/uj m for Fi = (a) (no phase 
noise), T 1 /2tt = 0.1 KHz (b), and r z /27r = 1 KHz (c). 
The cavity bandwidth has been fixed at ft/u; m = 0.5 and 
we have also fixed the center of the phase noise spectrum 
Q/2tt = 50 KHz. We see that laser phase noise has an 
appreciable effect on the log-negativity: for increasing Ti 
its maximum value decreases, and the parameter region 
where the steady state is entangled significantly narrows. 
Phase noise in particularly destructive close to the insta- 
bility threshold, where En is maximum when Ti = 0: as 
soon as Ti ^ 0, entanglement vanishes at the threshold 
for bist ability and the maximum value of En is achieved 
far from the threshold, still around A/uj m c± 1 and at 
intermediate values of the input power. In fact, in the 
presence of phase noise it is no more helpful to increase 
the input power because not only G but also \a s \ be- 
comes larger, amplifying in this way the effect of phase 
noise [see Eqs. (J^) and ( 10 i)]. 



Fig. 3 shows instead En versus the input power P 



and normalized cavity decay rate hz/uj m for Ti = (a), 
r^/27r = 0.1 KHz (b), and r z /27r = 1 KHz (c). The cav- 
ity detuning has been fixed at A/uj m = 1 and we have 
again fixed Q/2tt = 50 KHz. The destructive effects of 
laser phase noise manifest again by decreasing the maxi- 
mum achievable En and narrowing the parameter region 
with a nonzero entanglement. Again En vanishes at the 
bistability threshold, and for increasing phase noise en- 
tanglement is found only in the resolved sideband region, 

K/uJ m < 1. 

In Fig. 4 instead we study the dependence of entangle- 
ment upon the spectral properties of laser phase noise. It 
shows En versus the input power P and normalized cav- 
ity detuning A/uj m at fixed laser linewidth T^/27r = 0.1 
KHz, and for different values of the center of the fre- 
quency noise band, Q/2tt = 30 KHz (a), Q/2tt = 80 KHz 
(b), and Q/2tt = 140 KHz (c) (the bandwidth parame- 
ter 7 is correspondingly adjusted so that 7 = 0/2). The 
cavity decay rate has been fixed at K,/uj m = 0.5. The 
figures clearly show that the spectral properties of fre- 
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FIG. 4: (Color online) Contour plot of En versus the input power P and normalized cavity detuning A/uj m at fixed laser 
linewidth Fi/2tt — 0.1 KHz, and for different values of the center of the frequency noise band, Q/2tt — 30 KHz (a), Q/2tt — 80 
KHz (b), and £1/2tt = 140 KHz (c) (the bandwidth parameter 7 is correspondingly adjusted so that 7 = 0/2). The cavity 
decay rate has been fixed at n/uj m = 0.5. See text for the other parameter values. 
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FIG. 5: Contour plot of En versus the input power P and normalized cavity decay rate ft/cj m at fixed laser linewidth Ti/2ty = 0.1 
KHz, and for different values of the center of the frequency noise band, £}/2tt = 30 KHz (a), Q/2tt = 80 KHz (b), and Q/2tt = 140 
KHz (c) (the bandwidth parameter 7 is correspondingly adjusted so that 7 = Q/2). The detuning has been fixed at A/cj m = 1. 
See text for the other parameter values. 



quency noise have a strong effect on stationary entangle- 
ment, which in fact progressively worsens for broader and 
broader frequency noise spectrum. In fact entanglement 
is still considerable when Q/2tt = 30 KHz and 7/27r = 15 
KHz, but becomes extremely small when Vtj2i\ = 140 
KHz and 7/2-7T = 70 KHz, even with a moderate phase 
noise strength, Ti/2tt = 0.1 KHz. 

The relevance of the form of the noise spectrum is also 
evident in Fig. 5, where we show E^r versus the input 
power P and normalized cavity decay rate K,/cj m at fixed 
laser linewidth F//27T = 0.1 KHz, again for Q/2ir = 30 
KHz (a), n/2ir = 80 KHz (b), and ^/2tt = 140 KHz (c). 
The cavity detuning has been fixed at A/uj m = 1 and 
again the bandwidth parameter 7 is always correspond- 
ingly adjusted so that 7 = Q/2. Both the maximum 
achievable entanglement and the size of the parameter 
region with nonzero entanglement rapidly decreases with 
increasing bandwidth of the frequency noise spectrum. 



Again, entanglement is more robust against phase noise 
in the resolved sideband limit K/oJ m < 1. 



B. Approximate analytical expressions for En 

The exact expression of the logarithmic negativity 
stemming from the solution of Eq. (21) is cumbersome, 



but it is nonetheless possible to explain the above results 
by means of an approximate treatment which satisfacto- 
rily describes the effect of phase noise on En*. In fact, 
by following the approach of Ref. [29] , based on the solu- 
tion of Eqs. ([9| in the frequency domain, for each matrix 
element of the 4x4 stationary CM one can write 



SVijicu)}, 



(25) 



8 



where Vij(uj) refers to the situation with no laser noise 
and 5Vij (uj) is the correction due to the presence of laser 
phase noise. This latter correction is explicitly given by 



fiV iS (w) = 
where 

Xeff(w) = 



2H 2 | X effHI 2 ^MQHc i H* 

[k 2 + (lo- A) 2 ][k 2 + (w + A) 2 ] ■ 



(26) 



17mW ■ 



G 2 Aw n 



(k - iw) 2 + A 2 



(27) 

is the effective susceptibility of the mechanical oscillator 
modified by radiation pressure [50] , and c(uj) is the vector 



c(uj) = [GAu) m ,-iuGA,A(u)l 
(k - «cj)(cj^ - uj 2 - ij m uj)] T . 



As discussed in Ref. [50] (see also Ref. [53]), |x e ff (^)| 2 has 
a Lorenzian-like form, peaked at an effective frequency 
and with width given by an effective damping rate 



UJ. 



eff 



7^f , whose approximate expressions are given by [50] 



.eff 



G 2 Auu m (K 2 



•A 2 ) 



eff 



7n 



[k; 2 + (u; m - A) 2 ] [K? + (u; m + A)2] 
2G 2 Auj rn i<i 

[« 2 + (o; m - A)2] [«2 + (o; m + A)2] " 



(28 



(29) 



The modification of the mechanical frequency due to ra- 
diation pressure shown by Eq. (28) is the so-called "opti- 



cal spring effect" which may lead to significant frequency 
shifts in the case of low- frequency oscillators [54] . In the 
case of higher resonance frequencies (around 1 MHz) and 
red-detuned laser, as we are assuming here, the modifica- 
tion can be appreciable only very close to the bistability 
threshold. 

In a wide parameter region, and especially in the re- 
solved sideband regime n < uj mi the integrand in Eq. (26) 



is dominated by the peak of |x e ff(^)| 2 and therefore one 
can evaluate the corrections SVij by approximating the 
laser phase noise spectrum with a constant given by its 
value at the peak, S^(uj^). This means that in the re- 
solved sideband regime, the effect of laser phase noise 
on the stationary state of the system is completely de- 
scribed by a unique number, the phase noise spectrum at 
the effective mechanical frequency S^(uj^). Therefore 
the effects of laser phase noise can be minimized sim- 
ply by suppressing the noise spectrum at cj^f , and this 
fact explains why the explicit form of the phase noise 
spectrum is relevant. We have verified that the approx- 
imated CM obtained by replacing S\ (uj) with Si(lj^) 



in the phase noise terms of Eq. ( 26 ) satisfactorily repro- 
duces the behavior of En shown in Figs. 2-5, especially 
in the resolved sideband regime. 

Such an approximation for the laser phase noise con- 
tribution to the CM can also be used to derive an ap- 
proximate analytical expression for En in the parame- 
ter region very close to the bistability threshold. This 



regime is relevant for optomechanical entanglement be- 
cause, as first pointed out in Refs. [10] [29] and recently 
discussed in detail in Ref. [31 , En reaches its maximum 
value at the bistability threshold in the ideal case with 
no laser noise. As shown in Figs. 2-5, this is no more 
true in the presence of phase noise, when, on the con- 
trary, En always drops to zero close to threshold. This 
is well explained by the approximate expression of the 
symplectic eigenvalue rf~ which is obtained by approxi- 
mating S^(uj) with S^uj^ ) in Eq. (26), taking for G the 

threshold value G ~ y (k? + A 2 )cj m /A, and neglecting 
thermal noise terms. One gets 

1 U + bSi (cjeff) + cSi (cJeff) 2 + dSAuj^Y 



V2 1 



/ + #<S;(cJ eff ) 



(30) 



where 

a - 

b = 

C = 

d -- 

f ^ 

9 : 



«V + A 2 ) [4A 4 + 4A 2 {k 2 + u>l) + cui] ,(31) 
2\a s \ 2 A 2 K 2 [4(A 2 + /t 2 )(2A 2 + K 2 ) 
+6(A 2 + K 2 )u J 2 m +u J i], 
4K| 4 A 4 K [5(A 2 + / t 2 ) + 2u4], 
8k| 6 A 6 , 

8k 3 A V + A 2 ) (A 2 + k 2 + 5c4J , 



16|a s |VA 4 (A 2 + k 2 



(32) 
(33) 
(34) 
(35) 
(36) 



If laser phase noise is negligible, <S;(o; e ff) = 0, one gets 
(see also Ref. [3T] ) 
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4A 4 + 4A 2 (*; 2 +^ 2 n )+< 



16A 2 (A 2 + k 2 + 5ujI) ' 
This value is minimum at 

A/uj m 

at which En achieves its maximum value 



(37) 




E 



N 



In 




which can become at most In [5/3] ~ 0.51 in the resolved 
sideband limit n <C uj m [31] . 

In the presence of phase noise instead, the terms pro- 
portional to S^(uj e s) and its powers in Eq. (30) are always 
predominant because |a s | 2 is typically large, and there- 
fore r]~ may easily become larger than 1/2, which means 
no entanglement, even for not too large values of the 
laser linewidth Ti. This shows why entanglement van- 
ishes close to the bistability threshold as soon as phase 
noise is present, and also why the maximum En is at- 
tained at smaller input power for increasing phase noise. 
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C. Cooling 

The effect of laser phase noise on ground state cool- 
ing of the mechanical resonator has been already dis- 
cussed by various papers [33H36] . Ref. [35] in particu- 
lar provided an accurate estimation of the effect of laser 
noise and showed that the relevant quantity is just the 
frequency noise spectrum at the mechanical resonance, 
S^(uj m ): laser phase noise does not pose serious limita- 
tions to cooling provided that S^(uj m ) is not too large. 
The analysis of the effect of phase noise on the station- 
ary CM in the preceding subsection fully confirms such 
a prediction, i.e., we recover the results of Ref. [35] on 
ground state cooling, even though by means of a differ- 
ent treatment, based on QLE instead of the master equa- 
tion, and working in the frame rotating at the fluctuating 
frequency. 

In fact, the stationary mean energy of the mechanical 
oscillator is given by 



U 



■[(Sq 2 ) + (Sp 2 )} = nu m (n, 



eff 



(38) 



i.e., it is a linear combination of the stationary CM ma- 
trix elements Vu and V22 • Therefore one can again ex- 
ploit the fact that |Xeff(^)| 2 is strongly peaked, and ap- 
proximate the frequency integrals for (5q 2 ) and (Sp 2 ) by 
replacing S^u) with S^(oo^). Assuming also the weak 
coupling regime n ^> 7 m ?G and uj m ^> n7 m ,G which is 
the relevant one for ground state cooling, one obtains a 
result analogous to that of Ref. [35], which generalizes 
the expression of Refs. [50j [55j [56] by including the effect 
of laser phase noise, 



1 



7m -^op 



n7 n 



A , |^| 2 Ar oPc , /iCff> 



'SAO 



(39) 

Here we denote with r op = A- — A+ the net laser cooling 
rate, and 



In Fig. 6 we show n e ^ versus the input power P 
and normalized cavity detuning A/u m for Ti = (a), 
Ti/2n = 0.1 KHz (b), and Ti/2tt = 1 KHz (c). The cav- 
ity bandwidth has been fixed at ft/co> m = 1. Fig. 7 shows 
instead n e ff versus the input power P and normalized 
cavity decay rate hz/u m again for = (a), Ti/2tv = 0.1 
KHz (b), and Ti/2tt = 1 KHz (c). The cavity detuning 
is now fixed at A/uj m = 1, and in both figures we have 
fixed £}/2tt = 50 KHz. The effective occupancy is less 
affected than entanglement by laser phase noise: in fact 
for increasing Ti the parameter region where n e ff < 1 
becomes narrower but one can still reach ground state 
cooling for a realistic set of parameters, provided that the 
input power is not too large. Similarly to what happens 
for entanglement, ground state cooling is achievable only 
in the resolved side band limit for increasing laser noise 
strength Ti (see Fig. |7(c) ). 

In Figs. 8 and 9 instead we study the dependence 
of the phonon occupancy upon the spectral properties 
of laser phase noise. Fig. 8 shows n e ^ versus the in- 
put power P and normalized cavity detuning A/u m at 
fixed laser linewidth Ti/2tt = 0.1 KHz, fixed bandwidth 
ft/co> m = 1, and for different values of ^, Q/2tt = 30 
KHz (a), n/2ir = 80 KHz (b), and n/2ir = 140 KHz 
(c) (the bandwidth parameter 7 is always adjusted so 
that 7 = Cl/2). In Fig. 9 instead we plot n e ff versus the 
input power P and normalized cavity decay rate ft/c<j m , 
at fixed A/uj m = 1 and for the same set of values for 

and Q. As predicted by Eq. ([39]), it is just the noise 
spectrum at the effective mechanical resonance frequency 
which mainly affects cooling. When Q (and consequently 
7) is increased, the laser noise spectrum broadens and 
its value at cj^f increases as well. As a consequence, the 
parameter region where n e ff < 1 becomes narrower and 
narrower. In particular for a broader spectrum ground 
state cooling is better achieved at not too large values 
of input power and again in the resolved sideband limit 

K/uJ m < 1. 



kG 2 



2[^ + (A±u; m ) 2 ] 



(40) 



denotes the two scattering rates of laser photons into the 
Stokes sideband (A+) or the anti-Stokes sideband (A-) 
of the laser. Recall that, as pointed out in Ref. [50], 7 m + 
ff , i.e., it coincides with the effective mechanical 



op 



damping of Eq. (29). 



Eq. (|39j) reproduces Eq. (23) of Ref. [35] if we recall 
that in typical situations r op ^> 7 m and by restricting to 
the optimal condition for ground state cooling, i.e., the 
resolved sideband limit G, k <C u m and A « cj m , where 



the optical spring effect is negligible and ~ u m . We 



have compared the prediction of Eq. ( 39 ) with the exact 



solution of the stationary state of the system, and we 
have verified that it works very well in the relevant regime 
where ground state cooling is achievable. The behavior 
of n e ff for the same set of parameter values considered 
above is studied in Figs. 6-9. 



V. CONCLUSION 

We have studied the effects of the fluctuations of both 
the amplitude and the phase of the laser driving a cavity 
optomechanical system. We have analyzed the dynam- 
ics by adopting a quantum Langevin treatment, in which 
the phase noise dynamics has been included by means of 
additional auxiliary variables. We have linearized the dy- 
namics around the classical stationary state of the system 
and analyzed the dependence of the log-negativity of the 
stationary state and of the mechanical occupancy upon 
the various system parameters. We have also derived ap- 
proximate, but compact analytical expressions showing 
the effect of laser phase noise on these quantities. We 
have seen that, even though laser noise may have an ap- 
preciable affect on the quantum properties of the steady 
state, both cooling to the mechanical ground state and 
stationary optomechanical entanglement are still possible 
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FIG. 6: (Color online) Contour plot of n e s versus the input power P and normalized cavity detuning A/cj m for = (a), 
Ti/2ty = 0.1 KHz (b), and Yi/2n = 1 KHz (c). We have also fixed K/uj m = 1 and Q/2n = 50 KHz and the other parameter 
values are fixed in the text. 




FIG. 7: (Color online) Contour plot of n e ff versus the input power P and normalized cavity bandwidth ft/cj m for = (a), 
Ti/2tt = 0.1 KHz (b), and Ti/2tt = 1 KHz (c). We have also fixed A/cj m = 1 and Q/27T = 50 KHz and the other parameter 
values are fixed in the text. 
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FIG. 8: (Color online) Contour plot of n e ff versus the input power P and normalized cavity detuning A/o; m at fixed laser 
linewidth r^/27r = 0.1 KHz, fixed bandwidth n/uj m — 1, and for different values of the center of the noise spectrum, Q/2tt = 30 
KHz (a), Q/27T = 80 KHz (b), and Q/2n = 140 KHz (c) (the bandwidth parameter 7 is always adjusted so that 7 = £1/2). 
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FIG. 9: (Color online) Contour plot of n e ff versus the input power P and normalized cavity bandwidth A/u m at fixed laser 
linewidth Ti/2tv = 0.1 KHz, fixed detuning ft/cj m = 1, and for different values of the center of the noise spectrum, Q/2tv = 30 
KHz (a), Q/2-7T = 80 KHz (b), and Q/2tt = 140 KHz (c) (the bandwidth parameter 7 is always adjusted so that 7 = Q/2). 
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if state-of-the-art stable lasers are employed. (SQUALO project). 
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